Addendum to "Equation of state of classical Coulomb plasma mixtures" 
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Recently developed analytic approximation for the equation of state of fully ionized nonideal electron-ion 
plasma mixtures [Potekhin et al., Phys. Rev. E 79, 016411 (2009)], which covers the transition between the 
weak and strong Coulomb coupling regimes and reproduces numerical results obtained in the hypernetted chain 
(HNC) approximation, is modified in order to fit the small deviations from the linear mixing in the strong cou- 
pling regime, revealed by recent Monte Carlo simulations. In addition, a mixing rule is proposed for the regime 
of weak coupling, which generalizes post-Debye density corrections to the case of mixtures and numerically 
agrees with the HNC approximation in that regime. 

PACS numbers: 52.25.Kn, 05.70.Ce, 52.27.Gr 



I. INTRODUCTION 

The high accuracy of the linear mixing rule (LMR) for 
multicomponent strongly coupled Coulomb plasmas has been 
confirmed in a number of papers iH IH [H |j, HI |^ |7[l ■ Nev- 
ertheless, the accuracy of modern Monte Carlo (MC) calcu- 
lations allows one to reveal certain deviations from the LMR 
for the Coulomb energy U of binary ionic mixtures (BIM). 
On the other hand, for weakly coupled plasmas the Debye- 
Hiickel (hereafter DH) formula is applicable instead of the 
LMR. Several terms in the density expansion of U beyond the 
DH approximation were obtained by Abe and by Cohen 
& Murphy IgD (hereafter ACM) in the one-component plasma 
(OCP) case. 

In Ref. Jit], deviations from the LMR for BIM were stud- 
ied in the hypernetted chain (HNC) approximation and fitted 
by Pade approximants. In Ref. 1|4J, the LMR was confirmed 
by HNC method for polarizable background of partially de- 
generate electrons. In Ref. deviations from the LMR for 
strongly coupled BIM were studied using both HNC and MC 
techniques. The connections to the LMR for U were found 
to be of the same order of magnitude for HNC and MC, but 
numerically different; in particular, it does not depend on the 
mean ion Coulomb coupling parameter F according to HNC 
results, but decreases as function of F in MC simulations. 
These results were confirmed in Ref. where an analytic 
fit to the calculated corrections was suggested. The fitting for- 
mulae of Refs. L3i,2J are applicable only at F > 1; in particular 
they do not reproduce the DH Umit at F ^ (besides, the fit 
parameters in ^ are given only for 5 fixed ionic charge ratios 
from 2 to 8). 

In Ref. Ill], HNC calculations of BIM and three- 
component ionic mixtures (TIM) were performed in a wide 
range of values of F, charge ratios, and partial densities of the 
ion components, and a parametric formula was suggested to 
fit \ht fractional differences between the LMR and calculated 



plasma energies at any F in liquid multicomponent plasmas. 
It recovers the DH formula at F ^ 1 and gives a vanishing 
fractional difference from the LMR at F > L 

However, in the regime of strong coupling, the accuracy 
of the HNC method (typically a few parts in 1000, for U) is 
not sufficient to reproduce the values of the energies of mix- 
tures at the precision level needed to study deviations from 
the LMR (see, e.g., [5]). Indeed, according to Refs. |I1H0], 
these deviations are typically of the order of a few x(10~^- 
10^^) kT per ion (where k is the Boltzmann constant), while 
U FfcT per ion at F > 1. 

In this brief report, we suggest two improvements for ana- 
lytic treatment of ion mixtures. First, we introduce a mixing 
rule for weakly coupled plasmas, which provides an extension 
of the ACM formula to the case of ion mixtures and agrees 
with HNC results up to the values of the Coulomb coupling 
parameter F « 0.1 (whereas the DH approximation becomes 
inaccurate at F > 0.01). Second, using MC simulations of 
strongly coupled liquid BIM, supplementary to those already 
published in fsl S 01. we suggest a modified version of the 
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formula lllOll . which maintains the accuracy of the previous fit 
at intermediate and weak coupling, but delivers consistency 
with the MC data for strongly coupled Coulomb liquids. 

In Sec. we introduce basic notations and formulae; in 
Sec. Hn] we propose a mixing rule applicable at weak cou- 
pling; in Sec.|IV]we present a fitting formula for the internal 
energy of mixtures, applicable in the entire domain of F val- 
ues for weakly and strongly coupled classical Coulomb gases 
and liquids; and in Sec. [V] we summarize the results. 



n. BASIC EQUATIONS 

Let rie be the electron number density and rij the number 
density of ion species with charge numbers Zj (j=l,2,...). 
The total number density of ions is riions = X]j j ■ The elec- 
tric neutrality implies rig = (.^)nions- Here and hereafter the 
angular brackets denote averaging with statistical weights pro- 
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FIG. 1 : (Color online) Correction to the LMR Au = u — ulm as a 
function of T for BIM with ^2/^1 ^ 2, X2 = 0.2. HNC (crosses) 
and MC (dots) data are compared to the DH approximation (short- 
dashed lines), the modified ACM approximation (|5} (long-dashed 
lines), the fit from LI 0,1 (dot-dashed lines), and the present fit ^ 
(solid lines). 



portional to rij : 



{Z)=^x,Z. 



where xj = 



rii 



(1) 



The strength of the Coulomb interaction of ion species j is 
characterized by the Coulomb coupling parameter, defined (in 
CGS units) as Tj = {Zjef/ajkT = T^zf^, where = 
tteZy^ is the ion sphere radius, Fe = /a^kT, and Oe = 
(47r7ie/3)~^/'^. In other words, partial coupling parameters 
Vj and ion sphere radii aj are defined to be those of the OCP 
of ions of the jth kind at the same electron density Ue as in the 
considered multicomponent plasma. The Coulomb coupling 
in the mixture of different ions is conventionally characterized 
by the average coupling parameter T = Tf,{Z^/'^). 

A common approximation for the Coulomb contribution to 
the internal energy of a strongly coupled ion mixture is the 
LMR: 



1): 



(2) 



where u = U / NionskT is the reduced Coulomb energy, iVjons 
is the total number of all ions, and the subscript "LM" de- 
notes the linear-mixing approximation. Obviously, the LMR 
has the same form for the Coulomb contribution to the reduced 
Coulomb free energy / = F / NionskT . 

When the Coulomb interaction is sufficiently weak com- 
pared to the thermal energy, then the DH approximation 



FIG. 2: (Color online) Au = u — ulm as a function of F for BIM 
with Z2/Z1 = 2, X2 = 0.05. Here crosses (HNC^) correspond to 
All obtained from the HNC data using the OCP fit from llj for 
calculation of ulm, and asterisks (HNC^) correspond to A« from 
Ref. 1 5], where both u and mlm are based on the HNC results. Dots 
(MC^) correspond to Au calculated from the recent MC data for 
u, and itLM calculated from the OCP fit UliL while circles (MC^) 
represent MC data [5] for Aw. 



can be applied: mdh = {(f') /kTro-, where {q^) is the 
mean squared charge of the considered mixture and td ~ 
(/cT/47rnions('?^))^^^ is the Debye radius. For the model 
of ions in the "rigid" electron background, applicable if the 
electrons are extremely strongly degenerate, {q^) — e^{Z'^), 
whereas in the case of completely nondegenerate electrons, 
using our definition ([T]i of averaging over the ion species 
and taking into account the neutrality condition, we have 
{q^)=e\{Z^) + {Z)). 

In this paper we consider the model of rigid electron back- 
ground, but extension to the case of compressible background 
is possible by adjusting the parameter 5 in Eq. ( fTOl i below, ac- 
cording to the expression for {q'^). In Ref. [10 1 this extension 
was shown to be compatible with numerical HNC data for 
ion mixtures with allowance for electron polarization. 



III. WEAKLY COUPLED ION MIXTURES 

For a OCP at F < 1, a cluster expansion yields dll] 
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-F^/2(1.6875 \/3 InF - 0.23511) 



(3) 



where Ce = 0.57721 ... is the Euler constant. Here, the first 
term is the DH energy. 
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FIG. 3 : (Color online) The same as in Fig.Qbut for Z2/Z1 = 5 and 
X2 = 0.05. 



FIG. 4: (Color online) The same as in Fig.Qbut for Z2/Z1 = 8 and 
two values of X2'- 0.01 and 0.1 (marked near the dots). 



In order to generalize this expression to the case of multi- 
component Coulomb plasmas, let us write the OCP energy in 
the form 



u{T) = Tu{a/rD), 



(4) 



where a is the ion sphere radius for the OCP, and u is the 
Coulomb energy per ion in units of {eZY /a (u = —0.9 in the 
ion sphere model ifl^l ). Then the following relation holds in 
the DH approximation for multicomponent plasmas: 



{Z') ^ 



1/3 



(5) 



Let us assume that relation (|5]l can be applied also to the 
higher-order corrections beyond DH. In this case, according 
to Eqs. ([3]l and (|4|i, in the ACM approximation 



u{k) 



■ In K- 0.0149085 



ln(K) - 0.07369 



(6) 



Since k cx a/T for a fixed composition, / can be obtained 
from u by integration, which yields 



(7) 



where 



/» = -| - ^ (In 0.2263) 



-k' (0.02778 InK- 0.01946). 



(8) 



In Figs. [T]-@] deviations from the LMR, Au = u — ulm^ 
calculated according to Eqs. Q and (|6]l, are plotted by long- 
dashed lines and compared to the DH formula (short-dashed 
lines) and the HNC data (crosses). We see that suggested ap- 
proximation (|5]l agrees with the data to much higher F values 
than the DH approximation. 



IV. COULOMB LIQUIDS AT ARBITRARY COUPLING 



In order to find an analytic approximation for the correction 
to the LMR in the largest possible interval of F for ion gases 
and liquids, we have selected from the numerical HNC data 
IToll the subset related to F < 1, which counts 161 different 
combinations of Xi, Z2, and F in BIM and 54 combinations 
of X2, Z2, x-i, Zs, and F in TIM (assuming Zi = 1), sup- 
plemented this HNC data by numerical MC data for BIM at 
F > 1 (94 combinations of X2, Z2, and F), and looked for an 
analytic formula which provides a reasonable compromise be- 
tween simplicity and accuracy for representing this data. The 
MC data has been partly taken from the previous work [5,^,^ 
and partly obtained by new MC simulations using the same 
computer code as before. Our fitting formula for the addition 
to the reduced free energy / — F/NionskT, relative to the 
LMR prediction /lm, reads 



vl'^iz^'i^) 



LM 



V3 (l + aF")(l + 6F")/3 



(9) 
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where 6 is determined by the difference between the LMR and 
DH formula at T (exactly as in Ref. ifioll ): 



(5 = 1- 



(Z2)3/2 



for rigid electron background model, and 



5 = 



1)3/2) (^^2^ ^ ^^^)3/2 



(Z)l/2 (Z5/2 



(10a) 



(10b) 



for polarizable background. The expression jlObl ) for (5 is ex- 
act in the limit of nondegenerate electrons, but its use in Eq. 
(|9|l provides a satisfactory agreement with numerical data [4] 
obtained with allowance for the polarizability of partially de- 
generate electron gas (see [10]). 

The fit parameters a, b, and a are chosen so as to minimize 
the mean-square difference between the fit and the data for 
Au I ULM at r < 1 and for Au at F > 1, while the power index 
/3 is defined so as to quench the increase of A/ at F oo. 
These parameters depend on plasma composition as follows: 



2.6,5+1453 



(Z) 



2/5 



1 



6 = 0.0117 



(Z2) 



Z2)l/5' 

3 



/3 = 



2a 



(11) 



(12) 



The numerical difference of Eq. (|9]l from the formula in 
Ref. is small at F < 1, but at F > 1 the correction to the 
LMR prediction for the reduced internal energy 



Aw = F 



d{Af) 

dr 



3 aaF" 6a/3F" 
2 ^ 1 + aF" ~ 1 + 6F" 



A/ (13) 



now decreases at large F in agreement with the MC results. 
Moreover, Eq. (T3[ describes most of the data with much 
higher accuracy than the fit to Au suggested in Ref. Jltl for 
BIM at F > 1. 



A comparison of the numerical HNC data for A/ and Au 
and MC data for Au to Eq. (|9]l and to the previous fit ifToll 
shows that the present fit has nearly the same accuracy as the 
previous one for BIM at F < 1 (slightly worse for small 
Au/u, slightly better for larger Au/u), but it is generally 
better for TIM at F < 1 and substantially better for BIM 
at F > 1. Examples of F-dependences of Au are shown in 
Figs.[T]-3] where the dot-dashed lines correspond to the older 
fit and the solid lines to the present fit. The modification of the 
fit at small F values proves to be negligible, which has been 
checked by comparison of fractional differences between the 
Coulomb part of the free energy and the LMR prediction, as 
in Ref. |10], whereas the modification at large F can be sig- 
nificant, as confirmed by Figs. [T]-H] 



V. CONCLUSIONS 

We have reconsidered free and internal energies of classi- 
cal ionic mixtures in the liquid state, taking into account the 
results of HNC calculations in the regime of weak and moder- 
ate Coulomb coupling and MC simulations at strong coupling, 
and proposed two new analytic approximations for such mix- 
tures: the mixing rule (|5]l, which works well at the Coulomb 
coupling parameter F < 1, and the analytic fitting formula 
(|9]l, which is along with its derivative (fTSl l applicable at any 
values of F. 
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